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Abstract. - To probe the connection between the dynamic phase transition and stochastic 
resonance, we study the mean-field kinetic Ising model and the two-dimensional Josephson- 
junction array in the presence of appropriate oscillating magnetic fields. Observed in both 
systems are double stochastic resonance peaks, one below and the other above the dynamic 
transition temperature, the appearance of which is argued to be a generic property of the 
system with a continuous dynamic phase transition. In particular, the frequency matching 
condition around the dynamic phase transition between the external drive frequency and the 
internal characteristic frequency of the system is identified as the origin of such double peaks. 



When a system with an energetic activation barrier is weakly coupled to a temporally 
periodic driving, the inherent thermal stochastic noise can enhance the signal out of the 
system rather than weaken it. This behavior, called stochastic resonance (SR) |lj, has been 
studied extensively in systems with a few degrees of freedom; the appearance of a resonance 
peak at a finite temperature T$r, characteristic of SR, is conveniently explained by matching 
the two time scales: the deterministic one due to the external driving and the stochastic one 
which is inversely proportional to the thermal activation rate (Kramer rate) Recently, 
there have also been attempts to investigate the SR in more complex systems [0-^|. Here 
one should note that such a system with many degrees of freedom may undergo a phase 
transition at a nonzero temperature and have the characteristic time scale not described by 
the simple Kramer rate Q . Whereas the Kramer rate gives a diverging time scale only at zero 
temperature, the system with a finite-temperature continuous transition possesses a diverging 
characteristic time scale at a nonzero critical temperature T c , yielding finite values of the 
characteristic time both below and above T c . 

In this Letter we propose to extend the time-scale matching condition in such a way that 
the matching between the time scale of the external periodic driving and the characteristic 
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time scale associated with the phase transition determines the position of the SR peak. Since 
the time scales then match at two distinct temperatures, one below and the other above T c , it 
follows that the system with a continuous phase transition should show double SR peaks, one 
below and the other above the transition. In order to elucidate this scenario, we investigate 
the interplay between the dynamic phase transition and the SR behavior in the infinite-ranged 
kinetic Ising model and in the two-dimensional (2D) fully frustrated Josephson-junction array 
(FFJJA), driven by appropriate oscillating magnetic fields. Unlike the latter system, where 
the SR behavior has not been examined, the kinetic Ising model in an oscillating magnetic 
field has been studied previously via the mean-field approximation and numerical simulations, 
with regard to both the dynamic phase transition [^[j9| and the SR behavior |||4] , and double 
SR peaks have been found and discussed by Leung and Neda in Q|. However, the physics 
behind the emergence of double peaks has not been addressed per se: To our knowledge, 
neither the connection between the origin of the double peaks and the existence of a dynamic 
phase transition nor the generality of this connection has been recognized and addressed at 
all. 

We begin with the infinite-ranged kinetic Ising model described by the Hamiltonian 

j N 
H = --J2a i a j -h(t)J2^ (!) 

i<j i=l 

where the Ising spin (Tj at site i can have values ±1, J is the coupling strength, and the 
external oscillating magnetic field is given by h(t) = /i cos fit. Henceforth we set J = 1 and 
measure the temperature T in units of J/fcs. For the above infinite-ranged Hamiltonian, 
the mean-field approximation is exact and yields the equation of motion under the Glauber 
dynamics 

^=-m( f )+tanh( m(t) + ^), (2) 

where the magnetization m(t) is given by the ensemble average of the spin at time t. 

We first recapitulate the dynamic phase transition present in the system described by 
Eq. (||). In the absence of the external field, the system undergoes a thermodynamic phase 
transition at the transition temperature Tq = 1 between the low-temperature ferromagnetic 
phase (mo 7^ 0) and the high-temperature paramagnetic phase (mo = 0), where Too is the 
spontaneous magnetization in equilibrium. In the presence of the external oscillating field, on 
the other hand, the system displays a dynamic phase transition [^,[7] |oJ, for which the dynamic 
order parameter Q can be defined to be 

Q = lim Q n (3) 

71 — >OC 

with Q n denoting the time average during the nth period of external driving: 



^ /.27r(n+l)/n 

Q„ = — dtm(t). (4) 



2tt 



Note that Q separates the ferromagnetic phase (Q 7^ 0) from the paramagnetic phase (Q = 0) 
and determines the dynamic phase transition temperature T c . The phase boundary, obtained 
numerically by solving Eqs. (|J)-(|]) JsJ and analytically via the adiabatic approximation valid 
for sufficiently small f2 |lC| ], in general connects the points (T, h ) — (1,0) and (0, 1) on the 
T-h plane. For given frequency O, as the amplitude /i is increased from zero, the transition 
not only occurs at a lower temperature, lowering T c , but also changes its nature from a 
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Fig. 1 - Inverse of the relaxation time versus temperature, computed numerically from the behavior 
of Q n for /io =0.1 and Q,/2n = 0.1 (full line) and analytically via the linear response theory (dashed 
line). The relaxation time is shown to diverge, as the transition temperature is approached either 
from above or from below. 



continuous type to a discontinuous one. The division point at which nature of the transition 
changes moves along the phase boundary from (1, 0) towards (0, 1) as f2 is increased from zero. 
Good overall agreement exists in a broad range of temperatures between numerical results and 
analytical ones via the adiabatic approximation and the perturbation expansion JTo{ ] . 

A key quantity in our argument for the double SR peaks is the relaxation time r or 
cquivalcntly, the intrinsic frequency scale 1/r of the system, which describes the decay of the 
average magnetization: Q n — Q cx e~*/ r . Thus one may start from, e.g., the ordered state 
m(t=0) = 1, and extract the exponential decay of Q n to obtain numerically the relaxation 
time ||. At weak driving fields, the relaxation time can also be calculated analytically to the 
linear order in the amplitude h Q p|,|Io|]: 



1 _ 1 1 

t~ T cosh 2 (m /T)' 



(5) 



where the spontaneous magnetization mo is determined by the equation of state: mo = 
tanh(mo/T). Figure |l| displays the relaxation time for ho = 0.1 and = 0.1, determined 

by these two methods. It is shown that at the dynamical phase transition the character- 
istic frequency 1/r vanishes or equivalently, the relaxation time r diverges. Note also the 
good agreement between the analytical and numerical results, implying that Eq. (JsJ) gives an 
excellent approximation, except for the discrepancy in the temperature at which r diverges: 
Equation (^) leads to the divergence at the equilibrium transition temperature Tq = 1, whereas 
the actual relaxation time diverges at the dynamical transition temperature T c ss 0.98, slightly 
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Fig. 2 - Occupancy ratio versus temperature in the kinetic Ising model for ho — 0.1 and Q/2n = 0.1. 
Data directly from Eq. (^) together with numerical integration of Eq. (Q) are represented by small 
circles whereas the analytical expression given by Eq. ((?]) with r obtained numerically in Fig. [l] is 
plotted by the solid line. The dashed line corresponds to the linear response calculation [Eq. (fjj) with 
r given by Eq. (jj)]. Inset: the time scale matching condition with the solid line corresponding to the 
inverse relaxation time obtained numerically in Fig. ^ and the horizontal dashed line to the right-hand 
side in Eq. (^). The two crossing points below and above the transition temperature determine the 
positions of the two SR peaks in the main body. 

lowered by the external driving field. It is important here that a continuous dynamical tran- 
sition as well as a continuous equilibrium transition is associated with a divergent relaxation 
time. In particular the relaxation time can assume the same (finite) value at two temperatures, 
one above and the other below the transition. 

It is convenient to study the SR behavior in the system by means of the occupancy ratio 
(OR), defined to be the average fraction of the spins in the direction of the external field ||: 



At sufficiently low temperatures most spins are stuck to one of two spin directions (±1) and 
do not follow the direction preferred by the field h(t). Since hit) oscillates with period 2-it/Q,, 
most spins point in the same direction as h(t) during one half of the period while almost no 
spin follows h(t) during the other half; this results in OR « 1/2. At very high temperatures, on 
the other hand, each spin can take the values 1 and —1 with equal probability, again resulting 
in the value of OR close to 1/2, since about half the spins point in the same direction as h(t) 
at a given instant. Somewhere in between the spins may follow closely the external driving 
field, giving rise to a peak in the OR as a manifestation of SR. Within the linear response, 
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the OR takes the analytic form [nOL hi 



OR = - 
2 



— d> dt m(t) sgnfcos ilt) 
4?r / 



= l + h T - 1 (7) 

2 7T l + 2 r 2 ' K ' 

where § represents the integration over a full period of the driving field in the stationary 
state and r is given by Eq. (||). The resulting behavior of the OR with the temperature 
is shown in Fig. ||, which reveals most clearly the double-peak structure of SR. Figure ||, 
where the numerical result directly from the equation of motion (ph has also been displayed, 
demonstrates that the analytical result in Eq. ([?]) is rather accurate. Further shown in Fig. || 
is that the analytical result can even be improved by adopting in Eq. (j?]) the (more precise) 
numerical values of t shown in Fig. |l] instead of those given by Eq. (^|) . 

From Eq. (]?]), it is straightforward to obtain the peak position Tsr of the OR, which is 
determined by the condition (d/dT)OR = 0. This leads directly to the matching condition 
between the external and internal (intrinsic) time scales: 



r = 1 + Vl + sF*, (8) 

which reduces in the limit Q — > to the direct matching condition between the frequency 
scales: 

T^^n. (9) 

It is thus concluded that the resonance becomes strong when the external time scale I/O 
matches the intrinsic time scale r according to the matching condition in Eq. (||). Since r, 
diverging at T c , decreases away from T c , there in general exist two temperatures, one above 
and the other below the transition, at which Eq. (||) is satisfied for given fl: obviously, this 
gives rise to the double peaks. The inset of Fig. || displays the matching condition in Eq. (|J), 
the right-hand side of which is, for fl/2ir = 0.1, depicted by the horizontal dashed line: The 
two intersections determine precisely the positions of the two SR peaks in the main body 
of Fig. U, demonstrating explicitly the proposed matching condition for the existence of the 
double SR peaks. 

To disclose ubiquity of the double-peak structure associate with a dynamic phase transi- 
tion, we next consider the 2D L x L FFJJA, driven by a weak magnetic field, staggered in 
space and periodic in time, in addition to the uniform field introducing full frustration [see 
Fig. |^(a)]. Such a system is described by the Hamiltonian 

h = -jJ2 cos (<^ " h ~ ( 10 ) 

where J is the Josephson coupling strength, (pi denotes the phase of the superconducting order 
parameter at site i, and the summation is over all nearest neighboring pairs. The bond angle is 
given by the line integral of the vector potential and consists of two parts: Ay = A®j +SAij(t), 
where A^- describes the uniform magnetic field corresponding to half a flux quantum per 
plaquette (/ = 1/2) and JAy(t) takes into account the additional oscillating magnetic field 
of frequency O. When the oscillating field is absent (SA^j = 0), the ground state possesses 
the two-fold degeneracy and vortices form a 2x2 superlattice at zero temperature [ fL2"| . On 
the other hand, the presence of the oscillating field, which leads to the oscillating flux 8f(t) = 
/oCOsOi, lifts the two- fold degeneracy. For £1 -C 1, the energy splitting AE between the 
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Fig. 3 - (a) Schematic diagram of the fully-frustrated Josephson-junction array in the presence of 
an oscillating staggered magnetic field: /i,2 = 1/2 ± /o cos fit. Depending on the sign of cos fit, 
vortices energetically favor to be located on either shaded plaquettes or empty ones at given time; 
the double degeneracy present for fo = is lifted and the system can effectively be regarded as a 
two-state system in a double- well energy landscape as schematically shown in (b), where M denotes 
the staggered vorticity and AE is given in text, (c) A fully frustrated Josephson-junction array put 
on the top of a periodic array of magnetic particles, with either array driven periodically, provides a 
possible experimental realization. 



Fig. 4 - Occupancy ratio versus temperature in a fully frustrated two-dimensional Josephson-junction 
array driven by a weak oscillating staggered magnetic field. Simulation data for the 16 x 16 array at 
fo — 0.02 and = 0.01 are represented by squares, whereas the line is merely a guide to the eye. 
The double-peak structure is obvious although the peak below the transition is weaker than the one 
above. 



two states [each of which has the staggered vorticity ±1; see Fig. ||(b)] is given by pc[ : 
AE = 2 J-v/2 sin(7r/o cos fit/2). It is thus clear that this FFJJA subject to the oscillating 
staggered magnetic field contains all the basic ingredients required for SR behavior. Figure||(c) 
suggests one possible realization of this system in experiment. Periodic arrays of magnetic 
particles are already available and used frequently in experiment |L3| , and we believe that the 
relevant SR behavior can be observed through the use of various vortex imaging techniques. 

The dynamics of the system is described by the resistively-shunted junction model, stating 
the current conservation at each site: 

[<f>i - <f>j - Aij + sinO,; - - Aij) + Cij] = 0, (11) 

3 

where the prime restricts the summation to the neighbors of site i, and £y(i) is the random 
noise current across the junction with zero mean and correlation (Cij(t)Ckl(t')) = 2T(SikSji — 
8u8jk)5(t—t'). Note that energy (and temperature) has been measured in units of J and time 
in units of fi/Ae 2 JR with the shunt resistance R. 

To obtain the dynamic behavior, we integrate directly the set of equations of motion ( 111] ) 
and measure the staggered vorticity at each plaquette. The OR of the staggered vorticity, 
measuring how many plaquettes have the staggered vorticity preferred by the external driving, 
is then computed according to Eq. (||), with the oscillating flux f(t) taking the role of h(t) 
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Figure |4| displays the obtained OR as a function of the temperature T in the array of size 
L = 16 for fo — 0.02 and = 0.01. Observed clearly is the emergence of double SR peaks, one 
below and the other above the (dynamic) transition temperature T c w 0.443 p0[ , although 
the peak below the transition appears much weaker than the one above. Note that sufficient 
decrease of r is essential for the time-scale matching condition to be satisfied. In some systems, 
r may become large at low temperatures as well as near T c and the SR peak below T c may 
not be observed ]l4j]. 

In summary, we have proposed the extended time-scale matching condition for the SR, be- 
tween the external periodic driving and the relaxation time associated with a phase transition, 
thus elucidating the link between the appearance of double SR peaks and the existence of a 
dynamic phase transition. This has been demonstrated in the 2D fully frustrated Josephson- 
junction array as well as in the mean-field kinetic Ising model. The argument presented here 
is rather general, and we expect that the double-peak structure should appear ubiquitously in 
systems undergoing dynamic phase transitions with diverging relaxation time at the transition. 
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